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Abstract 

The phases of quantum interference effects in charged lepton production by neutrinos (neu- 
trino oscillations) following pion decay at rest and in flight, muon decay and nuclear j3— decay 
at rest as well as for 'muon oscillations' following pion decay at rest and in flight are calculated. 
■ The same phase is found for neutrino and muon oscillations following pion decay at rest. The 

results found disagree with the conventionally used value: 0i2 = Am 2 L/(2P). Differences to 
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previous treatments are briefly commented on. 
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This letter sketches briefly a derivation of spatially dependent interference effects (con- 
ventionally called 'neutrino oscillations') in the detection probability of charged leptons 
produced in neutrino interactions following the weak decays of unstable 'source' particles: 
pions, muons and (3— radioactive nuclei. Similar effects for the detection probability of 
the decay products of muons produced in pion decay ('muon oscillations') are also consid- 
ered. The calculations are based on Feynman's path amplitude formulation of quantum 
mechanics [1] in which the probability of transition from a set of initial states / = i\ 
to a set of final states F = Y^ m f™ is given by the relation: 
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where kj, j = l,n are (unobserved) intermediate quantum states. Only the essential 
elements of the calculation and the results are presented here. Full details are given 
elsewhere [2]. 

The application of Eqn(l) to pion decay at rest is illustrated in the ideal experiment 
shown in Fig.l. A 7r + comes to rest in a stopping target T at time to as recorded by the 
counter Ca (Fig. la)). The pion at rest constitutes the initial state of the path amplitudes. 
In Fig. lb) and Fig.lc) are shown two alternative histories of the stopped pion. In Fig. lb) 
the pion decays at time t± into the neutrino mass eigenstate \ui >, of mass mi, and in 
Fig.lc) into the neutrino mass eigenstate |z/2 >, of mass m2, at the later time £2- If 
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mi > m 2i then, for a suitable choice of the times t\ and t 2 , interference between the path 
amplitudes corresponding to the different physical processes shown in Fig. lb) and Fig.lc) 
will occur when a neutrino interaction (z/ 1; v 2 )n — > e~p takes place at time and distance 
L from the pion decay point (Fig. Id)). The final state of the neutrino interaction event is 
also that of the path amplitudes. For this experiment, the path amplitudes corresponding 
to the two alternative histories of the decaying pion are: 

At = < e~p\T R \riUi > U ei D(x f - x u t D - t h m^U^ 

< Vifi + \T R \7r + > e-^t^O^-to.m,) i = 1,2 (2) 

Here U ai is the Maki-Nakagawa-Sakata (MNS) [3] matrix element giving the charged 
current coupling strength of the charged lepton of flavour a (a — e, /i, r) to the neutrino 
mass eigenstate v { . The 'reduced' matrix elements < e~p\T R \nUi > and < ViH + \T R \it + > 
are given by setting U a i = 1 in the charged current, so that, for example, < e~p\T\nVi >= 
U e i < e~p\T R \nvi >. Since the purely kinematical effects of non- vanishing neutrino masses 
are expected to be negligible, then, to a very good approximation, \v\ > and \ u 2 > may be 
replaced, in the reduced matrix elements, by the massless neutrino wavefunction |z/ >. 
The former are then independent of neutrino flavour. In the case of a 2 x 2 MNS matrix 
for the first two generations of leptons, the unitarity of the matrix implies that it is 
completely defined by a single real parameter, conventionally chosen to be an angle, 9, 
such that: 

U el = U le = U„ 2 = U 2lx = cos9 (3) 
U e2 = U 2e = -U lll = -U lll = wy9 (4) 

In Eqn(2), m n and r„- are the pion pole mass and decay width and D is the Lorentz 
invariant configuration space propagator [4, 5] of a neutrino or the pion. In the limit of 
large time-like separations, or of on-shell particles, appropriate to the experiment shown in 
Fig.l, D ~ exp[— imAr] where m is the pole mass of the particle and At the increment 
of proper time corresponding to the path. In the following, the additional functional 
dependence ~ (to/ At) 2 of D in the asymptotic region (leading to solid angle correction 
factors) is neglected. With these approximations: 

D(Ax, At, to) ~ exp[-iTO V / (At) 2 - (Ax) 2 } 
= exp[— imAr] 

= exp[-iA0] (5) 

The phase increments, A<f>, corresponding to the paths of the neutrinos and the pion in 
the amplitudes are: 
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Atf = m t Ar t = -f At, = -±L (6) 

Acpi = m n (ti - t ) = m n {t D - t ) 

~ m n (t D - t ) - m n L ^1 + ^ j (7) 
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where TOj, E h P i and i>j are the mass, energy, momentum and velocity of the mass 
eigenstate ui and 
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The neutrinos are assumed to follow classical rectilinear trajectories such that At-i = 
L/vi = EiL/Pi and the time dilatation formula At = 7 At = EAr/m has been used 
in Eqn(6). In Eqn(7) the neutrino velocity" Vi is expressed in terms of the neutrino 
mass to order m 2 . Using Eqns(5-7), neglecting terms of 0(m 4 ) and higher, and replacing 
the massive neutrino wavefunctions in the reduced matrix elements by those of massless 
neutrinos, the path amplitudes of Eqn(2) may be written as: 



Ai = <e p\T R \nu >U ei U ill <u fi + \T R \7i + > 
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where the neutrino times-of- flight t{ 1 — to — U have been introduced and 

0o = m w (L — t D + to) 
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Using now Eqn(l) to calculate the transition probability, and integrating over the de- 
tection time to [2], gives, for the probability to observe the reaction (z/!,z/ 2 )n — > e~p at 
distance L from the decay point: 



P(e~p\L) = C N (v;ir) sin 2 9 cos 2 9(1 - F U (T„) cos^f) 



(11) 
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and Eqns(3) and (4) have been used. Here Cn(u;h) is an L independent normalisation 



factor and Am 2 



TOT 



■to 2 .. The first and second terms in the second member of Eqn(12) 



are the contributions of A0J 1 " and A0^ to the interference phase. The latter is a factor 
of two larger than in the conventional value [6, 7] (f)±2 = Am 2 L/(2P). A<pJ gives a 
numerically large (m 7T /2Po = 2.34) contribution to the oscillation phase. 

The above calculation is readily repeated for the case of detection of muon decay 
fi + — > e + (i/i, i^X^i, v<i) at distance L from the n + decay point. The phase increments 
analagous to Eqns(6) and (7) are: 
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"Units with h — c — 1 arc used. 
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The result found for the time-integrated decay probability is: 



P(e + uV\L) = C Ar (//;7r)(l-F^(r 7r ) sin 29 cos <ft£) 



(16) 



where 
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The first and second terms in the second member of Eqn(17) are the contributions of 
A(f)J^ and A(f)f to the interference phase. It is interesting to note that neutrino and 
muon oscillation phases are the same for given values of Am 2 and L. For oscillation phases 
<Pi2 = 012 = 1) ^e damping factors of the oscillation term, due to the non- vanishing 
pion lifetime take the values F v (Y n ) = l-1.58x lfr 16 and F>*(T n ) = l-4.4x 1(T 16 . This 
damping effect is thus completely negligible in typical neutrino oscillation experiments 
with oscillation phases of order unity. 

Formulae like Eqn(ll) have been derived in a similar manner [2] for '77 M — > 77 e ' os- 
cillations b following muon decay at rest, detected via the process (Vi,V 2 )p — > e + n and 
l V e — > Ve oscillations following nuclear (3— decay, detected via the same process. The 
results found, for the time-integrated decay probabilities, are: 



P(e + n, n\L) = Cn{v\ aO sin 2 9 cos 2 9 
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In these formulae, P v is the momentum of the detected V\,V2 , Ep is the total energy 
release in the (3— decay process, and damping corrections due to the finite lifetimes of the 
decaying particles are neglected. 

Finally, in Ref. [2], the cases of neutrino and muon oscillations following the decay in 
flight of ultrarelativistic pions were considered. The oscillation probability formulae are 
the same as Eqns(ll) and (16) respectively, with the phases: 



#f(in flight) 
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The neutrinos or muons are detected at distance L from the decay point, along the 
direction of flight of the parent pion, with decay angles 9 V and 9 ^ respectively. 

Because the parent pion and the daughter muon in the decay process n —> jiu are un- 
stable particles, their physical masses W n and have distributions depending, through 
Breit-Wigner amplitudes, on their pole masses m n and m M and decay widths and T^. 



b li the neutrinos are massive, the 'lepton flavour eigenstates' V e and v ^ are unphysical. However the phrase '77^ — > v e 
oscillations' is still a useful and compact way of describing experiments where anti-neutrinos, produced in association 
with a muon, give rise to a detection event containing an electron. The amplitudes for all physical processes contain, 
however, only the physical states 
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Energy-momentum conservation in the pion decay process then leads to a distribution of 
path amplitudes with different neutrino momenta, an effect neglected in the above dis- 
cussion. In Ref.[2] corrections resulting from the (coherent) momentum smearing due to 

and (incoherent) smearing due to W n are calculated using a Gaussian approximation 
for the Breit-Wigner amplitudes. The resulting damping corrections due to W n to the 
interference terms in Eqns(ll) and (16) are found to be vanishingly small. For Am 2 = 
(leV) 2 and L = 30m (typical of the LNSD [8] or KARMEN [9] experiments) the corre- 
sponding damping factors are found to be I — 1.3 x 10~ 29 (for neutrino oscillations) and 
1 — 5.2 x 10~ 30 (for muon oscillations) [2]. The damping corrections from the variation 
of Wfj, are even smaller. 

Corrections to Eqns(ll) and (16) due to thermal motion of the decaying pion and finite 
target and detector sizes have also been evaluated in Ref.[2]. For neutrino oscillations 
with Am 2 = (leV) 2 and L = 30m and a room-temperature target, the damping factor of 
the interference term is found to be 1 — 6.7 x 10~ 10 and the shift in the oscillation phase 
to be 1.2 x 10~ 9 rad. In summary, all known sources of damping of the interference terms 
in Eqn(10) and (15) are expected to be completely negligible in any forseeable neutrino 
or muon oscillation experiment. 

The first published calculation of the neutrino oscillation phase [10] gave a result (only 
the contribution of the neutrino propagators was taken into account) in agreement with 
Eqn(6) above. A later calculation [11] assumed instead that the phase of the neutrino 
propagator evolves as Et, i.e. according to the non-relativistic Schrodinger equation. 
The assumptions were also made that the different neutrino mass states are produced at 
the same time and have equal momenta. The first assumption leads to the 'standard' 
formula for the oscillation phase: </>i 2 = Am 2 L/(2P) [6, 7] which has subsequently been 
used for the analysis of all neutrino oscillation experiments. The Lorentz-invariant phase 
of Eqn(6) ~ m 2 t/E evidently agrees with the result of Ref.fll] in the non-relativistic 
limit where E ~ m, but such a limit is clearly inappropriate to describe experiments with 
ultra-relativistic neutrinos. 

The most important difference in the treatment given in the present paper to previous 
ones that have appeared in the literature is allowing the possibility for the different 
neutrino mass eigenstates to be produced at different times. Only in this way can the 
constraints, of both space-time geometry (the detection event is at a unique space-time 
point), and exact energy- momentum conservation in the decay process [12], be satisfied. 

The other new feature is the inclusion of the important contribution to the oscilla- 
tion phase from the propagator of the decaying particle, a necessary consequence of the 
different production times of the different mass eigenstates. 

In the standard treatment, the common production time of the mass eigenstates follows 
from the assumption that a 'neutrino flavour eigenstate', that is a coherent superposition 
of mass eigenstates, is produced in the decay of the source particle. It has been shown 
that the production of such states in pion decay is incompatible with the measured 
branching ratio r(7r — > ev)/T{it — > fiv) and the values of the MNS matrix elements 
deduced from atmospheric and solar neutrino oscillations [13]. Actually, as first pointed 
out by Shrock [14, 15], as a consequence of the Standard Model structure of the charged 
lepton current: 

j^CCyept = ^ q7m( i _ 75 )[4^ (23) 

a,i 
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the different neutrino mass eigenstates are produced incoherently in different physical 
processes. This is the basic assumption of the calculations presented above. Clearly, in 
this case, the different mass eigenstates may be produced at different times in the (clas- 
sically) alternative histories corresponding to the different path amplitudes. It has also 
been recently demonstrated [16] that the factor of two difference in the contribution of 
neutrino propagation to the oscillation phase between the calculation presented above and 
the standard formula is a necessary consequence of the assumption of equal production 
times for all mass eigenstates and hence equal space-time (as contrasted to kinematical c ) 
velocities. In the same paper [16] it is also shown that different kinematical hypotheses 
that lead to different, but unequal, kinematical velocities (equal momenta, equal ener- 
gies, or exact energy-momentum conservation in the production process) differ in their 
predictions for the oscillation phase only by terms of 0(mf), and so are all equivalent at 

OK). 

It may be remarked that the physical interpretation of 'neutrino oscillations' provided 
by the path amplitude description is different from the conventional one in terms of 'lep- 
ton flavour eigenstates'. In the latter, the amplitudes of different lepton flavours in the 
neutrino are supposed to vary harmonically as a function of time. In the amplitudes 
for the different physical processes in the path amplitude treatment there is, instead, no 
variation of lepton flavour in the propagating neutrinos. Only in the detection process 
itself, the associated MNS matrix elements give different coupling strengths between the 
different neutrino mass eigenstates and the final state charged lepton and the interfer- 
ence effect occurs between the different path amplitudes that is described as 'neutrino 
oscillations'. In the case of the observation of the recoil muons no such 'lepton flavour 
projection' occurs in the detection event, but exactly similar interference effects are pre- 
dicted to occur. As previously emphasised [17], the 'flavour oscillations' of neutrinos, 
neutral kaons and b- mesons are just special examples of the universal phenomenon of 
quantum mechanical superposition that is the physical basis of Eqn(l). 

The most important practical conclusion of the work presented here is that identical 
information on neutrino masses is given by the observation of either neutrinos or muons 
from pion decay at rest. In view of the possibility of detecting muons simply and with 
essentially 100 % efficiency, in contrast to the tiny observable event rates of neutrino 
interactions, the recent indications for 'z/ M — > 17 e ' oscillations with Am 2 ~(leV) 2 following 
fi + decay at rest [8] could be easily checked by a search for muon oscillations following 
7r + decay at rest. For Am 2 ~(leV) 2 , the first absolute maximum of the interference term 
in Eqn(16) occurs at L ~ 8m. Note that Eqn(16) is valid for any muon detection process 
that does not distinguish between events where v x or u 2 are produced. 

c The 'kinematical' velocity, Vki n is defined = p/E. 
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Figure 1: The space-time description of an experiment in which neutrinos produced in the processes 
7r + — > (J, + (vi, v 2 ) are detected at distance, L, via the processes (z/ 1; z/ 2 )n — > e~p (see text). 
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